(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     38818,       1062]
NotebookOptionsPosition[     35965,       1003]
NotebookOutlinePosition[     36393,       1020]
CellTagsIndexPosition[     36350,       1017]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"kernelU", "[", "u_", "]"}], ":=", 
  RowBox[{"1", "/", "2"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"kernelT", "[", "u_", "]"}], ":=", 
  RowBox[{"1", "-", 
   RowBox[{"Abs", "[", "u", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"GammaT1", "=", 
   RowBox[{"(", GridBox[{
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{"kernelT", "[", "u", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{"u", " ", 
          RowBox[{"kernelT", "[", "u", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]},
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{"u", " ", 
          RowBox[{"kernelT", "[", "u", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "^", "2"}], " ", 
          RowBox[{"kernelT", "[", "u", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]}
     }], ")"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"PsiT1", "=", 
   RowBox[{"(", GridBox[{
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"kernelT", "[", "u", "]"}], "^", "2"}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{"u", " ", 
          RowBox[{
           RowBox[{"kernelT", "[", "u", "]"}], "^", "2"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]},
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{"u", " ", 
          RowBox[{
           RowBox[{"kernelT", "[", "u", "]"}], "^", "2"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "^", "2"}], " ", 
          RowBox[{
           RowBox[{"kernelT", "[", "u", "]"}], "^", "2"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]}
     }], ")"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"varpsilonT1", "=", 
   RowBox[{"(", GridBox[{
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "^", "2"}], " ", 
          RowBox[{"kernelT", "[", "u", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]},
      {
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "^", "3"}], " ", 
          RowBox[{"kernelT", "[", "u", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"u", ",", "0", ",", "1"}], "}"}]}], "]"}]}
     }], ")"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.838981276432822*^9, 3.8389812969673524`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"610fa782-a012-40ac-9a3d-0b26ad929317"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP1", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", "9.809037"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", "6.924716"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.4105", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.4105", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.96891114"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "1199"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
  3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
  3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
  3.8389818979241486`*^9, 
  3.8389819053863187`*^9}},ExpressionUUID->"cb932476-559d-4879-9499-\
ba5148423968"],

Cell[BoxData["0.3344486518587857`"], "Output",
 CellChangeTimes->{{3.8389817836338205`*^9, 3.8389817952531047`*^9}, 
   3.8389819204597597`*^9},ExpressionUUID->"6f39d80e-2ce7-4d19-9775-\
13774e22567f"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP2", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", "96.22992"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "113.4154"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"6.1524", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"6.1524", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.50600542"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "2620"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
  3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
  3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
  3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
  3.838982038708338*^9}},ExpressionUUID->"c0fa7ada-de8c-4500-9347-\
d8fffe44b6bb"],

Cell[BoxData["0.17321069845931153`"], "Output",
 CellChangeTimes->{
  3.8389820480665317`*^9},ExpressionUUID->"5800899b-2006-47c6-84f4-\
5992dc426465"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP3", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "3.556184"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "6.357568"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{".4698", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{".4698", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.50585306"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "17304"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
  3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
  3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
  3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
  3.838982119634227*^9}, {3.838982463964139*^9, 3.838982466042178*^9}, {
  3.8836649899280157`*^9, 3.8836649900161686`*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"2f694cb9-d711-4150-a2ed-dde19da7b29d"],

Cell[BoxData["0.2384647652313094`"], "Output",
 CellChangeTimes->{
  3.8389821280702486`*^9, 3.83898246723377*^9, {3.8836649305652413`*^9, 
   3.883664958337183*^9}, 3.8837298630135016`*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"cc831170-d7f8-44bf-b6c5-47b3d5a42e2a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP4", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "1.288178"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "3.678939"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.9576", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.9576", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.51388908"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "22558"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 
   3.838982494665718*^9},ExpressionUUID->"a1bf4cfb-156f-456e-b571-\
a4802b0c0e01"],

Cell[BoxData["0.3193503834380909`"], "Output",
 CellChangeTimes->{
  3.838982498649441*^9},ExpressionUUID->"025555b2-d8be-494c-86e0-\
1d082dd90411"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP5", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "15.81248"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", "15.99251"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.5724", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.5724", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "1.0470558"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "1066"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 
   3.8389826332838993`*^9}},ExpressionUUID->"367f3324-5e63-4817-beba-\
2222233a6633"],

Cell[BoxData["0.14742752908187218`"], "Output",
 CellChangeTimes->{
  3.838982658999728*^9},ExpressionUUID->"f3d7acba-5397-4240-98fd-\
721c01fc1a85"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP6", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", "0.4000707"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "0.7679749"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.0788", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.0788", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.10597359"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "5161"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}},ExpressionUUID->"dee6efc7-de1e-4e9a-8dc6-\
9a874cf95c01"],

Cell[BoxData["0.2884579515067119`"], "Output",
 CellChangeTimes->{
  3.838982785593503*^9},ExpressionUUID->"42fb0451-57a4-450a-a792-\
87884511599a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP7", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "8.654692"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", "1.088098"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.544", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.544", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.57018627"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "12894"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}, {3.838982876097715*^9, 
   3.8389829321670046`*^9}},ExpressionUUID->"072b800c-9b2e-4dea-b37b-\
fbb1e063251f"],

Cell[BoxData["0.15904777111017804`"], "Output",
 CellChangeTimes->{
  3.838982934823831*^9},ExpressionUUID->"145788c2-8c29-414a-8238-\
9325d1f7c10d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP8", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "1.789442"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", "0.1467769"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.2167", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.2167", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.53387165"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "12527"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}, {3.838982876097715*^9, 3.8389829321670046`*^9}, {
   3.8389830326951714`*^9, 3.8389830329629254`*^9}, {3.8389830685687613`*^9, 
   3.838983107419041*^9}},ExpressionUUID->"50cb4fcc-da8b-4494-8204-\
768e26743aa6"],

Cell[BoxData["0.21406577795022147`"], "Output",
 CellChangeTimes->{
  3.8389831097326765`*^9},ExpressionUUID->"a822aab6-73ca-427e-85a4-\
788786ce8d38"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP9", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", 
     RowBox[{"-", "21.28845"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", "18.03547"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.7691", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.7691", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "1.8107637"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "4742"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}, {3.838982876097715*^9, 3.8389829321670046`*^9}, {
   3.8389830326951714`*^9, 3.8389830329629254`*^9}, {3.8389830685687613`*^9, 
   3.8389831729154015`*^9}, 
   3.8389832982426553`*^9},ExpressionUUID->"c5555f7b-c06b-480f-ae03-\
bc81da96da2f"],

Cell[BoxData["0.1013463204347533`"], "Output",
 CellChangeTimes->{
  3.8389833076080837`*^9},ExpressionUUID->"76034609-16c8-43a9-affe-\
588e7252e07a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP10", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", "3.12233"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "0.3569933"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{".0942", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{".0942", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.93443796"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "2627"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}, {3.838982876097715*^9, 3.8389829321670046`*^9}, {
   3.8389830326951714`*^9, 3.8389830329629254`*^9}, {3.8389830685687613`*^9, 
   3.8389831729154015`*^9}, 3.8389832982426553`*^9, {3.838983342974842*^9, 
   3.8389833436958933`*^9}, {3.8389833792710266`*^9, 3.838983413314373*^9}, 
   3.8389834791810217`*^9},ExpressionUUID->"69774cde-9eb8-42a5-84bd-\
facd89155b3d"],

Cell[BoxData["0.1482731746323293`"], "Output",
 CellChangeTimes->{3.8389834194542494`*^9, 
  3.8389834829017534`*^9},ExpressionUUID->"6048324a-4a28-4e31-abb7-\
8d89d5b54e00"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "DGP11", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"beta2Yl", "=", "2.586736"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"beta2Yr", "=", 
     RowBox[{"-", "0.2259673"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yr", "=", 
     RowBox[{"0.4556", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sigma2Yl", "=", 
     RowBox[{"0.4556", "^", "2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"f0", "=", "0.64426993"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"NObs", "=", "39406"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BT1", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"beta2Yr", "-", "beta2Yl"}], ")"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "varpsilonT1"}], 
         ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"VT1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"sigma2Yr", "+", "sigma2Yl"}], ")"}], "/", "f0"}], " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Inverse", "[", "GammaT1", "]"}], ".", "PsiT1", ".", 
          RowBox[{"Inverse", "[", "GammaT1", "]"}]}], ")"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ConsT1", "=", 
     FractionBox["VT1", 
      RowBox[{"4", " ", 
       RowBox[{"BT1", "^", "2"}]}]]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"hT1", "[", "n_", "]"}], ":=", 
    RowBox[{
     RowBox[{"ConsT1", "^", 
      RowBox[{"(", 
       FractionBox["1", "5"], ")"}]}], 
     RowBox[{"n", "^", 
      RowBox[{"(", 
       RowBox[{"-", 
        FractionBox["1", "5"]}], ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{"hT1", "[", "NObs", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.838981411115714*^9, 3.8389814459338427`*^9}, {
   3.8389815195369453`*^9, 3.8389815987922707`*^9}, {3.8389816632259226`*^9, 
   3.838981735296054*^9}, {3.8389817938949327`*^9, 3.8389817944489512`*^9}, {
   3.8389818979241486`*^9, 3.8389819306668806`*^9}, {3.8389820043593884`*^9, 
   3.838982119634227*^9}, {3.838982349234395*^9, 3.8389823493050747`*^9}, {
   3.838982403899104*^9, 3.838982451713699*^9}, 3.838982494665718*^9, {
   3.83898258941508*^9, 3.8389826332838993`*^9}, {3.8389826952567625`*^9, 
   3.8389827665779924`*^9}, {3.838982876097715*^9, 3.8389829321670046`*^9}, {
   3.8389830326951714`*^9, 3.8389830329629254`*^9}, {3.8389830685687613`*^9, 
   3.8389831729154015`*^9}, 3.8389832982426553`*^9, {3.838983342974842*^9, 
   3.8389833436958933`*^9}, {3.8389833792710266`*^9, 3.838983413314373*^9}, {
   3.8389834486425686`*^9, 3.838983504586606*^9}, {3.838983729231306*^9, 
   3.8389837586831303`*^9}},ExpressionUUID->"09de0cbd-e062-4b6e-a323-\
af9c72d911e0"],

Cell[BoxData["0.19006411266045528`"], "Output",
 CellChangeTimes->{
  3.8389837705184574`*^9},ExpressionUUID->"b16dfcee-8a8e-4744-9577-\
3685871c8c4d"]
}, Open  ]]
},
WindowSize->{1587, 968},
WindowMargins->{{0, Automatic}, {Automatic, 134}},
Magnification:>2.1 Inherited,
FrontEndVersion->"12.3 for Microsoft Windows (64-bit) (July 9, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"2a3ef0f8-eca7-4a40-b186-f0ca65d512dd"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 3148, 93, 355, "Input",ExpressionUUID->"610fa782-a012-40ac-9a3d-0b26ad929317"],
Cell[CellGroupData[{
Cell[3731, 117, 2396, 66, 581, "Input",ExpressionUUID->"cb932476-559d-4879-9499-ba5148423968"],
Cell[6130, 185, 201, 3, 69, "Output",ExpressionUUID->"6f39d80e-2ce7-4d19-9775-13774e22567f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6368, 193, 2465, 67, 581, "Input",ExpressionUUID->"c0fa7ada-de8c-4500-9347-d8fffe44b6bb"],
Cell[8836, 262, 151, 3, 69, "Output",ExpressionUUID->"5800899b-2006-47c6-84f4-5992dc426465"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9024, 270, 2605, 69, 581, "Input",ExpressionUUID->"2f694cb9-d711-4150-a2ed-dde19da7b29d"],
Cell[11632, 341, 269, 4, 105, "Output",ExpressionUUID->"cc831170-d7f8-44bf-b6c5-47b3d5a42e2a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11938, 350, 2615, 70, 581, "Input",ExpressionUUID->"a1bf4cfb-156f-456e-b571-a4802b0c0e01"],
Cell[14556, 422, 148, 3, 69, "Output",ExpressionUUID->"025555b2-d8be-494c-86e0-1d082dd90411"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14741, 430, 2643, 70, 581, "Input",ExpressionUUID->"367f3324-5e63-4817-beba-2222233a6633"],
Cell[17387, 502, 149, 3, 91, "Output",ExpressionUUID->"f3d7acba-5397-4240-98fd-721c01fc1a85"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17573, 510, 2696, 70, 589, "Input",ExpressionUUID->"dee6efc7-de1e-4e9a-8dc6-9a874cf95c01"],
Cell[20272, 582, 148, 3, 91, "Output",ExpressionUUID->"42fb0451-57a4-450a-a792-87884511599a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20457, 590, 2745, 71, 616, "Input",ExpressionUUID->"072b800c-9b2e-4dea-b37b-fbb1e063251f"],
Cell[23205, 663, 149, 3, 91, "Output",ExpressionUUID->"145788c2-8c29-414a-8238-9325d1f7c10d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23391, 671, 2850, 72, 616, "Input",ExpressionUUID->"50cb4fcc-da8b-4494-8204-768e26743aa6"],
Cell[26244, 745, 151, 3, 91, "Output",ExpressionUUID->"a822aab6-73ca-427e-85a4-788786ce8d38"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26432, 753, 2877, 73, 616, "Input",ExpressionUUID->"c5555f7b-c06b-480f-ae03-bc81da96da2f"],
Cell[29312, 828, 150, 3, 91, "Output",ExpressionUUID->"76034609-16c8-43a9-affe-588e7252e07a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29499, 836, 3001, 74, 616, "Input",ExpressionUUID->"69774cde-9eb8-42a5-84bd-facd89155b3d"],
Cell[32503, 912, 174, 3, 91, "Output",ExpressionUUID->"6048324a-4a28-4e31-abb7-8d89d5b54e00"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32714, 920, 3081, 75, 616, "Input",ExpressionUUID->"09de0cbd-e062-4b6e-a323-af9c72d911e0"],
Cell[35798, 997, 151, 3, 124, "Output",ExpressionUUID->"b16dfcee-8a8e-4744-9577-3685871c8c4d"]
}, Open  ]]
}
]
*)

